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O ' Entanglement of formation for a class of higher dimensional quantum 

mixed states is studied in terms of a generalized formula of concurrence 
Oh. for A^-dimensional quantum systems. As applications, the entangle- 

i ment of formation for a class of 16 x 16 density matrices are calculated. 
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Quantum entanglement plays important roles in quantum communication, information 
processing and quantum computing [1], such as in the investigation of quantum teleporta- 
tion [2, 3, 4], dense coding [5], decoherence in quantum computers [1] and the evaluation of 
quantum cryptographic schemes [6]. To quantify entanglement, a well justified and mathe- 
matically tractable measure of entanglement is needed. A number of entanglement measures 
such as the entanglement of formation and distillation [7, 8, 9], negativity [10, 11], and rel- 
ative entropy [9, 12] have been proposed for bipartite states [6,8,12-17]. Nevertheless, most 
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proposed measures of entanglement involve extremizations which are difficult to handle 
analytically. 

The "entanglement of formation" is intended to quantify the amount of quantum commu- 
nication required to create a given state [7]. Although it is defined for arbitrary dimensions, 
so far no explicit analytic formulae for entanglement of formation have been found for sys- 
tems larger than a pair of qubits, due to the fact that two dimensional bipartite mixed 
states are special in many ways [16], except for some special symmetric states [17]. 

In this letter we study the entanglement of formation for a class of higher dimensional 
quantum mixed states. For certain A'"-dimensional pure quantum systems, we show that the 
entanglement of formation is a monotonically increasing function of a kind of generalized 
concurrence. As applications, the entanglement of formation for a class of 16 x 16 density 
matrices is calculated in detail. The method applies to a large class of quantum states. 
The construction of these states are presented for N dimensional, N — 2^^+^, 2 < k E IN, 
bipartite systems. 

Let H. be an A^-dimensional complex Hilbert space with orthonormal basis Cj, i = 1, A^. 
A pure state on 7i is generally of the form, 

N 

with normalization 

N 

E = 1 ■ (2) 

The entanglement of formation E is defined as the entropy of either of the two sub-Hilbert 
space H®n [8], 

= -TV (pi log2 pi) = -TV (p2 log2 P2) , (3) 
where pi (resp. P2) is the partial trace of over the ffist (resp. second) Hilbert space 

Let A denote the matrix with entries given by in (1). pi can be expressed as 

Pi = AAl (4) 

For a given density matrix of a pair of quantum systems onTi^Ti, consider all possible pure- 
state decompositions of p, i.e., all ensembles of state of the form (1) with probabihties 
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p^^Pi\A){Al ^Pi^^ 

i=l i=l 

for some M G IV. The entanglement of formation for the mixed state p is defined as the 
average entanglement of the pure states of the decomposition, minimized over all possible 
decompositions of p, 

M 

E{p)^ mm J2PiEm)- (5) 

i=l 

It is a challenge to calculate (5) for general N. Till now a general explicit formula of 
E{p) is obtained only for the case N — 2. In this case (3) can be written as 



where 

h{x) = -x\og2X- (1 -a;)log2(l -x), 
C is called concurrence [15]: 

C(IV')) = = 2|aiia22 - 012021!, 

where lijj) = ay (g) ay\ip*), is the complex conjugate of 1-0), ay is the Pauli matrix, 



ay 



As £^ is a monotonically increasing function of C, C can be also taken as a kind of measure 
of entanglement. Calculating (5) is reduced to calculate the corresponding minimum of 
C(p) = min J2iLiPiC{\ipi)), which simplifies the problems. 

For N >3, there is no such concurrence C in general. The concurrences discussed in [18] 
can be only used to judge whether a pure state is separable (or maximally entangled) or not 
[19, 20]. The entanglement of formation is no longer a monotonically increasing function 
of these concurrences. Nevertheless, for a special class of quantum states, we can find 
certain quantities (generalized concurrence) to simplify the calculation of the corresponding 
entanglement of formation. 

[Theorem 1]. If AA^ has only two non-zero eigenvalues (each of which may be degenerate), 
the maximal non-zero diagonal determinant D of AA^ is a generalized concurrence. The 
entanglement of formation of the corresponding pure state is a monotonically increasing 
function of D. 



[Proof]. Let Ai (resp. A2) be the two non-zero eigenvalues of AA^ with degeneracy n 
(resp. m), n + m < N. That is, 

D = A^A^ . (6) 
From the normahzation of \iJj), one has Tr{AA^) = 1, i.e., 

nAi + 171X2 = 1 • (7) 

Ai (resp. A2) takes values (0, ^) (resp. (0, ^)). In this case the entanglement of formation 
of 1'^) is given by 

= logs Ai - mA2 logs A2 . (8) 

According to (6) and (7) we get 

dD ~ 1-nXi- mAi m J mAi ' ^ ^ 

which is positive for Ai e (0, ^). Therefore E{\t/j)) is a monotonically increasing function of 
D. is a generalized concurrence and can be taken as a kind of measure of entanglement 
in this case. ■ 

Remark: We have assumed that Ai, A2 7^ in our theorem. In fact the right hand 
side of (9) keeps positive even when Ai (or equivalently A2) goes to zero. Hence E(\t/j)) 
is a monotonically increasing function of D for Ai G [0, ^] (resp. A2 G [0, ^]) satisfying 
the relation (7). Nevertheless if Ai = (or A2 = 0), from (6) one gets D = 0, which 
does not necessarily mean that the corresponding state 1-0) is separable. As £'(|'0)) is just 
a monotonically increasing function oi D, D only characterizes the relative degree of the 
entanglement among the class of these states. 

From (7) and (8), the quantum states with the measure of entanglement characterized 
by D are generally entangled. They are separated when n = 1, Ai — > 1 (A2 — > 0) or m = 1, 
A2 ^ 1 (Ai 0). For the case n — m> 1, all the pure states in this class are non-separable. 
In this case, we have 

E{\ij)) = n f-a;log2a;- {- - x)\og2{- - x)] , (10) 
V n n J 



where 



1/1 1 



and 



d = 2nD^ = 2n^AiA2. (11) 

In this case wc define d to be the generahzed concurrence, d takes value from to 1. From 
(10) one can show that E{d) is a convex function (that is, curving upward): 

Instead of calculating E{p) directly, one may calculate the minimum decomposition of 
D{p) or d{p) to simplify the calculations. In the following, as an example, we calculate the 
entanglement of formation for a class of mixed states with N — A. 

We consider a class of pure states (1) with the matrix A given by 



/ b ai bi \ 

—b Ci di 

Qi Ci — e 

\ bi di e J 



(12) 



ai,bi,ci,di,b,e G C. The matrix AA^ has two eigenvalues with degeneracy two, i.e., 
n = m = 2. 

det{AA^) = |6ici - aidi + be\^. (13) 
According to our theorem, the generalized concurrence 

d ^ A\biCi - aidi + be\ (14) 
is a kind of measure of entanglement for all pure states of the form (12). 

Let p be a 16 X 16 matrix with only non-zero entries pi^ie = P2,i5 = ~P3,u — P4,io — 

P5,12 = P6,ll = P7,13 = -P8,8 = -P9,9 = PlO,4 = Pll,6 = Pl2,5 = Pl3,7 = -Pl4,3 = Pl5,2 = Pw,l = 

1. d in (14) can be written as 

d^mpr)\ = \mm, m 

where {{ilj\ip)) — {'^\p'4>*). 

We now calculate the entanglement of formation for a special class of mixed states. Let 
^ denote the set of pure states (1) with A given as the form of (12). We consider all mixed 
states with density matrix p such that its decompositions are of the form 

M M 

P = T.P^\^^){'^^l E^^ = 1' 1^0 ^ (16) 
1=1 1=1 



Let s < 16 be the rank of p and \vi), i = l,...,s, be a complete set of orthogonal 
eigenvectors corresponding to the nonzero eigenvalues of p, such that {vi\vi) is equal to 
the ith eigenvalue. Other decomposition {ju^j)} of p can then be obtained through unitary 

transformations: 

K)=Et^^K), (17) 

where U is a. t x t unitary matrix, t > s. The states \wi) are so normalized that p = 
Yli \wi){u'i\- It is obvious that for any G \E', the state of complex linear combination of 
iV'j) (unitary transformations) also belongs to ^. 

The decomposition according to the orthogonal eigenvectors \vi) of p is not the one 
satisfying (5) in general. As the generalized concurrence can be written as the form (15) for 
a pure state, we consider the quantity {{wi\wj)). From (17) we have 

where the matrix r is defined by Tij = {{vi\vj)). The matrix p in (15) is a symmetric 
one. Therefore r is also symmetric and can always be diagonalized by a unitary matrix U 
such that UtU'^ — diag{Ai, A^) [21]. The diagonal elements Aj, in deceasing order, can 
always be made to be real and non-negative. Since Utt*W is also diagonal, Aj are just the 
square roots of the eigenvalues of rr*. It is straightforward to check that they are also the 
eigenvalues of the Hermitian matrix R = ^ \fpVP*V\fP-, or, alternatively, the square roots of 
the eigenvalues of the non-Hermitian matrix ppp*p. 

Hence there always exits a decomposition consisting of states i = 1, . . . ,s, such 
that 

{{wi\wj)) = KSij. (18) 
We can now deal with the problem in a way similar to [15]. Set 

|yi) = ki), \yj)^i\wj) for j = 2, (19) 
Any decomposition can be written in terms of the states \yi) via the equation 

\zi) = i:v*j\yj), 

where V is at x s matrix whose s columns are othonormal vectors. 
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The average concurrence of a general decomposition is given by 

(d) = Y.iivyy'^u = l:\l:iv^J)%J > (20) 



where Y is the real diagonal matrix defined by Y^j = {{yi\yj))- Using the fact that 
\{Vijf\ = 1, one gets 

16 

{d)>\T.(Vij)%j\>^i-E^i- 

ij i=2 

Therefore the minimum decomposition of the generalized concurrence is 

16 

d{p)^A,-J2K. (21) 

i=2 

Similar to the case N — 2, there are decompositions such that the generahzed concurrence 
of each individual state is equal to d{p). Therefore the average entanglement is E{d{p)). 

Different from the case N = 2, the entanglement of formation of density matrices (16) 
can not be zero in general. As every individual pure state in the decompositions is generally 
an entangled one, this class of mixed states are not separable. 

In the following we call an A^-dimensional pure state (1) d-computable if A satisfies the 
following relations: 

det{AA^) = {[Aurr/^ 

22 

det{AA^ - A/djv) = (A^ - ||yl||A + [A][A]*)^/'', 

where [A] and \\A\\ are any quadratic forms of (these quadratic forms could be different 
for different matrix A), Idj^ is the N x N identity matrix. We denote A the set of matrices 
satisfying (22), which implies that for A & A, AA^ has at most two different eigenvalues 
and each one has order N/2. Formula (21) can be generalized to general A^^ x N'^ density 
matrices with decompositions on A^-dimensional d-computable pure states. 

A class of A^-dimensional, N = 2'',2 < k E IN, ^-computable states has been constructed 
in [22] . These states give rise to a special class of density matrices with decompositions in 
these pure states, and the entanglement of formation for these density matrices can be 
calculated analytically according to the method above. 

Let yl be an A" X A" matrix with entries Ojj e C, i,j = 1,...,N, with the following 
properties: 



Set 

a —c\ 

c d ' 
/ 

where a,c,dEL C For any 61, ci e C, a 4 x 4 matrix A4 E A can be constructed in the 
following way, 



B2 A2 

— 4* r** 



f bi a -c \ 

—bi c d 

—a —c — Ci 

c —d ci y 



where 



B2 = 61J2, C2 = Ci J2, -^2 



1 
-1 



satisfies the relations in (22): 

^441 = [(6,ci + ad + c2)(6ici + ad + c2)*]2 = ([^4][>l4]*)2, 

A44 - Xldi = (A^ - (6161 + ciCi + aa* + 2cc* + drf*)A 
+ (61 ci + ad + c^) (61 ci + ad + c^)*)^ 

= (A2-P4||A+[A4][A4]*)^ 

where 

[^4] = (61C1 + ad + c^), 11^411 = 6161 + CiCi + aa* + 2cc* + dd*. 



& A can be obtained from A4, 



/ S4 A4 \ 



—^4 04 



where 



-B4 — 62<^4) C4 — C2J4, J4 



62, C2 e C 



(23) 



(24) 



(25) 



(26) 



/ 1 \ 
10 
0-100 

\-i y 

For general construction of high dimensional matrices ^2^+1 ^ A, 2 < k E IN, we have 

/ A2k \ I &fct-^2'= ^2*= 



\ (-1)^4, y 1^ {-l)'^Al, cJl, 



(27) 



J2k + 1 



( J2u\ 



:28i 



where 6fc, Cfc e C, ^2*; = hkJik, — CkJ2k. It can be verified that ^2* satisfies the following 

relations [22]: 

\A2k,.Al^,\ = {[A2k+.][A2k,.rf = [((-l)^6,c, - [A2kM-l)'^blcl - 

1^2*+! ^2-^+1 - A/d2fc+l| = {^^ - \ \A2k+l\\X+ [A2k+l][A2k+l]*y'' . 

(29) 

Therefore the states given by (27) are d-computable. In terms of (11) the generahzed 
concurrence for these states is given by 

d2k+i = 2'=+^|[^fe+i]| = 2'=+>feCfe + 6fe_iCfc_i + ... + bici + ad + c^]. 

Let P2k+i be a symmetric anti-diagonal 2^^^^ x 2^'^"'"^ matrix with all the anti-diagonal 
elements 1 except for those at rows 2*^+^ -1 + 5(2*^+^ - 2), 2*^+^ + s(2'=+2 - 2), 2*^+2 - 1 
s(2'=+2 _ 2), 2'=+2 + s(2'=+2 - 2), s = 0, 2^=+^ - 1, which are -1. d2k+i can be written as 

d2k+l = |(^2fe+l|P2fe+l^2fe+l)l = l((V^2fe+l|^2fc+l))|, (30) 

where 

2k+l 

lV'2*+i) = Yl {A2k+i)ij ei (g) Cj. (31) 

According to the calculations on entanglement of formation for d-computable states, for a 
22fc+2 ^ 22'=+2 density matrix p2k+2 with decompositions on pure states of the form (31), its 
entanglement of formation is given by E{d2k+i{p22k+2)), where 

(^2^+1 (p22*+2 ) — ^1 ~ ^i, (32) 

i=2 

and Qj, in decreasing order, are the the square roots of the eigenvalues of the matrix 

P22k+2p2k+l P22fe+2P2fe+l ■ 

We have studied the entanglement of formation for a class of higher dimensional quan- 
tum mixed states. It is shown that for certain A'"-dimensional pure quantum systems, the 
entanglement of formation is a monotonically increasing function of a generalized concur- 
rence. From this generalized concurrence the entanglement of formation for a large class of 
quantum states can be calculated analytically. The physical properties of these states are 
remained to be studied further. 
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